1 



(N 
O 
O 
(N 

o 

Q 



Recent Progress in Perturbative Quantum Field Theory* 

Zvi Bern'*''^ 

'^Department of Physics and Astronomy 
University of California at Los Angeles 
Los Angeles, CA 90095 

In this talk, the recent breakthrough in two-loop perturbative calculations is reviewed, with emphasis on the 
applications to phenomenological studies. The recent precision measurement of the anomalous magnetic moment 
of the muon by the Brookhaven g — 2 collaboration is also used as an illustrative example of the importance of 
Standard Model precision perturbative calculations to further our understanding of Nature. 
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1. INTRODUCTION 

Precision calculations in the Standard Model 
have a long history dating back to Schwinger's 
original QED calculations. The importance of 
precision Standard Model calculations in modern 
day particle physics was recognized by the 1999 
Nobel Prize awarded to Gerard 't Hooft and Mar- 
tinus Vcltman. The tools and theoretical founda- 
tion provided by 't Hooft and Veltman have led 
to a plethora of precision calculations, including 
the famous 'blueband' plot, shown many times at 
this conference, bounding the mass of a Standard 
Model Higgs boson, mn :S 196-230 GeV at 95% 
CL§. 

In this talk we describe recent progress in pre- 
cision calculations, as well as present examples 
illustrating the importance of higher loop compu- 
tations for comparing theory to experiments. The 
most spectacular recent example of this is the pre- 
cision measurement of the anomalous magnetic 
moment of the muon by the Brookhaven g — 2 col- 
laboration This experiment is an ideal 
example due to its impressive sensitivity to the 
quantum corrections arising from all three com- 
ponents of the Standard Model: Quantum Elec- 
trodynamics, Electroweak Theory, and Quantum 
Chromodynamics. The incredible agreement of 
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approximately seven digits of accuracy between 
theory and experiment is a testament to the 
remarkable experimental and theoretical efforts 
that have gone into this attempt to discover new 
physics beyond the Standard Model. 

This talk also focuses on the very recent break- 
through in our ability to calculate two-loop quan- 
tum corrections when more than a single kine- 
matic invariant is present. This breakthrough 
is the culmination of many years of theoretical 
effort. To illustrate the promise of this break- 
through, we present two very recent applications 
improving our understanding of Higgs physics, 
which is of central importance in particle physics. 
As one example, just two weeks before the confer- 
ence, Anastasiou and Melnikov |^ computed the 
exact NNLO result for next-to-next-to-leading or- 
der (NNLO) QCD corrections to inclusive Higgs 
production, in perfect agreement with the ear- 
lier series expansion result of Harlander and Kil- 
gore 0. The second example, completed a month 
prior to the conference, is giuon fusion into a pho- 
ton pair jsj^. This is a significant component 
of the background to the preferred search mode 
for discovering and measuring the Higgs mass at 
the Large Hadron Collider, if the Higgs is light 
{AIh < 140 GeV), as suggested by precision elec- 
troweak measurements. 

These examples represent initial steps in ap- 
plying the breakthrough to phenomenology. It 
is now clear that over the next few years, NNLO 
calculations of jet observables at both hadron and 
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lepton colliders will become available ||10|,y_l|jl2| . 
(Other related issues may be found in Frix- 
ione's QCD theory talk.) Precision cal- 
culations of Bhabha scattering is another ex- 
ample that will surely benefit from the break- 
through (l|,|l|. This process is important for 
measuring luminosity at e"'"e~ colliders. More 
generally, many observables which could not pre- 
viously be calculated to the required order of per- 
turbation theory necessary to match experimen- 
tal precisions will now be computed. 

This talk is organized as follows. In Section ^ 
the theoretical calculations required for a preci- 
sion prediction of the anomalous magnetic mo- 
ment of the muon are reviewed briefly. In Sec- 
tion 1^ some examples illustrating the clever tricks 
that enter challenging perturbative calculations 
are given. A brief overview of the status of loop 
calculations is presented in Section ^. In Section ^ 
we discuss the recent breakthrough in our abil- 
ity to compute quantities of interest with more 
than a single kinematic variable. In Section ^ 
two very recent examples of applications of the 
new advances to phenomenology are presented. 
Prospects for the future are described in Sec- 
tion 0. 

2. EXAMPLE: ANOMALOUS MAG- 
NETIC MOMENT CALCULATIONS 

The anomalous magnetic moment of the muon 
is an especially interesting quantity because of its 
sensitivity to high scale new physics. It is about 
40,000 ~ (m^/TOe)^ times more sensitive than 
the electron magnetic moment in most potential 
scenarios of physics beyond the Standard Model. 
The incredible precision of both the theoretical 
predictions and the experiments, coupled with its 
sensitivity to potential new physics makes mea- 
surement of the anomalous magnetic moment of 
the muon an important means of searching for 
new physics beyond the Standard Model (see e.g. 
refs. [[60). 

The magnetic moment of an elementary parti- 
cle is related to it spin via the gyromagnetic ratio: 
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Figure 1. The Feynman diagram describing the 
Dirac prediction for the magnetic moment of the 
muon. 
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Figure 2. The Feynman diagram containing the 
leading quantum correction to the magnetic mo- 
ment of a lepton, first obtained by Schwinger. 



From the Dirac equation, the gyromagnetic ratio 
is g — 2 for spin 1/2 particles such as the muon 
or the electron. The tree-level Feynman diagram 
describing the Dirac prediction is depicted in Fig- 
ure ^ Quantum mechanical effects represented 
by loop diagrams such as the one in Figure ^, 
however, alter this result very slightly leading to 
a non-vanishing 'anomalous' magnetic moment: 
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Computations of anomalous magnetic moments 
have a long history dating back to Schwinger's 
original calculation of the 0{a) contribution in 
1948 (l|. 

The incredibly precise experimental result 
announced at this conference by Yannis Se- 
mcrtzidis on behalf of the Brookhaven muon 
g ~ 2 collaboration is 



al^P = 11 659 204(7) (5) x 10 



-10 



(3) 



which is a significant improvement over the pre- 
vious measurement [|], 



^cxp 



11659 202(15) X 10 



-10 



(4) 



where /2 is the magnetic moment and S the spin. 



and an even greater improvement on the series of 
experiments at CERN that ended in 1977. 
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A discrepancy between the experimental result 
and the theoretical prediction using the Standard 
Model can in principle indicate new physics be- 
yond the Standard Model, 



(5) 



cxp ^SM new physics 

In Figure |^ a sample diagram which could lead 
to a New Physics contribution to the anomalous 
magnetic moment is shown. The virtual u parti- 
cle depicted in the diagram would be some as yet 
undiscovered particle coming from a New Physics 
scenario, such as supersymmetry. 
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dard Model. We now briefly review the intense 
theoretical effort that has gone into producing 
Standard Model predictions with the requisite 
precision. More thorough reviews may be found 
elsewhere [p^ 19 5|. 



2.1. QED contributions 

Calculations of the QED contributions to the 
anomalous magnetic moment of the muon and 
other Icptons have a long history (see e.g. 
refs. ]20| , ^ , p2|) starting with the seminal work 
of Schwingcr |l8| . The QED effects are by far the 
largest of the contributions, but are under the 
best theoretical control. 



Figure 3. New physics can lead to a contribution 
to the anomalous magnetic moment beyond that 
found in the Standard Model. The dotted line 
labeled by v represents a particle of a New Physics 
scenario. 



The initial Brookhaven g result caused some 
excitement in the community because of a 
2.6(7 discrepancy between theory and experiment. 
However, as described below in the subsection on 
hadronic contributions, a correction was made in 
the theoretical computation that significantly re- 
duces the discrepancy. Amusingly, the improved 
experimental results announced at this conference 
bring back the discrepancy to essentially the orig- 
inal level. However, the theoretical uncertainties 
associated with the hadronic corrections are still 
the subject of intense debate and it is not yet 
clear if the discrepancy is due to new physics. 

For the purposes of this talk what is impor- 
tant is not the small discrepancy between theory 
and experiment, but the astonishing agreement to 
better than six significant digits. This agreement 
is a great triumph for modern particle physics, 
from both the experimental and theoretical sides. 

On the theoretical side the computation re- 
quires a firm grasp of all components of the Stan- 



Figure 4. Examples of two-loop QED contribu- 
tions to the anomalous magnetic moment. 



The combined results of many years of effort is 
that the QED contribution is 

OED ^ ^+ 0.765857376(27) f-V 

+24.05050898(44) (6) 

-M26.07(41)(^)'^ 

-1-930(170) (^)V--- 

These results are from the theory update of Czar- 
necki and Marciano |l^, compiling the results of 
refs. [p3|j2^ . The four-loop (a/Tr)^ contributions, 
for example, consist of a total of 891 Feynman 
diagrams which were evaluated by Kinoshita and 
his collaborators |2^ via numerical means. (A 
recent correction actually shifts this very slightly, 
but with negligible effect on the prediction of the 
muon anomalous magnetic moment p4| .) The 
five loop {a/n)^ corrections are not exact calcu- 
lations but are instead estimated using various 
approximations pq]. An important contribution 
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which significantly reduces the uncertainty is the 
exact analytical result for the three-loop compu- 
tation . The net result for the QED contri- 
bution to the anomalous magnetic moment of the 
muon is: 



^QED ^ 11658 470.59(29) x 10" 



(7) 



Although the QED contribution is the largest of 
the contributions, the theoretical uncertainty is 
tiny compared to the experimental uncertainty in 
eq. (H) due to the high order of these calculations. 

2.2. Electroweak contributions 

The electroweak contribution is much smaller 
than the QED contribution because of the sup- 
pression due to W or Z propagators. Neverthe- 
less, the Brookhaven experiment is sensitive to 
these tiny corrections. The one-loop electroweak 
corrections displayed in Figure ^ were computed 
about 30 years ago . The more recently com- 
puted two-loop corrections are surprisingly 
large, shifting the one-loop results by about 23%. 
This large shift is due mainly to the appearance 
of large logarithms of the form ln(Mf /M^). The 
two-loop electroweak contribution is non-trivial 
to calculate since it is described by a total of 
1650 Feynman diagram of which about 200 are 
important. Two sample diagrams are depicted 
in Figure |^. The final result for the electroweak 
contributions including the two loop corrections 



af^ = 15.1(4) X 10 



-10 



(8) 



The theoretical uncertainty in this quantity is 
tiny compared to the uncertainty in the exper- 
imental measurement However, the over- 
all contribution of the expression is well within 
the sensitivity of the experiment, providing an 
exquisitely indirect confirmation of electroweak 
theory |,|. 

2.3. Hadronic contributions 

The most problematic of the contributions are 
the hadronic corrections to the anomalous mag- 
netic moment. The underlying difhculty is that 
Quantum Chromodynamics is strongly interact- 
ing, making it much more difficult to carry out 
computations. At the low energies relevant for 
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Figure 5. One-loop electroweak Feynman dia- 
grams contributing to the anomalous magnetic 
moment of the muon. 



7 W 
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Figure 6. Two of the 1650 two-loop electroweak 
Feynman diagrams contributing to the anomalous 
magnetic moment of the muon. 



the muon anomalous magnetic moment, the QCD 
coupling becomes large and perturbation theory 
breaks down. As yet, there are no available meth- 
ods for performing a first principles computation 
of the hadronic contributions, though perhaps in 
the future lattice gauge theory may be able to 
provide some useful input. Currently, experimen- 
tal input and model calculations are used. 

The two basic categories of hadronic contribu- 
tions are the vacuum polarization and light-by- 
light contributions depicted in Figure |7| and |. 
The more significant vacuum polarization contri- 
butions shown in Figure |^ can be obtained from 
measurements of the annihilation cross section 
e~^e~ — > hadrons, by making use of the optical 
theorem and dispersion relations (see, for exam- 
ple, refs. iQ). The light-by-light contributions 
depicted in Figure ||, however, can only be ob- 
tained via model calculations with relatively large 
theoretical uncertainties (29|. However, since the 
latter contributions are much smaller, the vacuum 
polarization uncertainty dominates. The over- 
all value of the light-by-light contributions, has 
recently been dramatically modified because of 
the discovery of a sign error in previous calcu- 
lations [p0|; this correction helps reduce the dis- 
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crepancy between theory and experiment. 

The status of the uncertainties in the hadronic 
corrections was reviewed at this conference by 
Teubner §. Since the time of the conference a 
number of new theoretical analyses have also ap- 
peared ]3l|| . Depending on various assumptions, 
currently the discrepancy is somewhere between 
1.6 and 3(7. It will, however, be some time before 
there is a final consensus on its significance. 




Figure 7. The hadronic vacuum polarization con- 
tribution. The wavy lines are photons and the 
blob represents hadronic contributions. 



3.1. Helicity Methods 

As a first example, consider the fivc-giuon tree 
amplitude described by Feynman diagrams of 
which two are depicted in Figure |. These dia- 
grams are straightforward to evaluate by direct 
algebraic evaluation on a computer, using one 
of the algebraic programming languages such as 
FORM @, MAPLE, and so forth. The result of 
evaluating the Feynman diagrams is shown in Fig- 
ure |lO|, using a microscopic font in order to illus- 
trate its apparent complexity. The result shown 
in this figure is actually just the coefficient of one 
of the color factors; the complete answer is much 
larger. 



o 




Figure 9. Two of the Feynman diagrams describ- 
ing the five-gluon tree amplitude. 



Figure 8. The hadronic light-by- light contribu- 
tion. 



3. CHALLENGING PERTURBATIVE 
COMPUTATIONS 

The area of higher order perturbative computa- 
tions in quantum field theory is a rather challeng- 
ing field, requiring clever ideas and algorithms for 
dealing with the difficulties. Instead of present- 
ing a systematic description of the methods used 
in state-of-the-art computations, here we provide 
a selection of examples to illustrate the type of 
ideas that can be helpful in challenging pertur- 
bative calculations. More complete descriptions 
may be found in various review articles, such as 
refs. ||,||. 



After some careful thought, a number of au- 
thors concluded that the mess of the type shown 
in Figure |l^ is not due to any inherent property 
of gluon scattering but is only due to using an in- 
efficient representation of the amplitude | ]35| , ^ . 
A much better representation is to use helicity 
states instead of polarization vectors. For each 
color a gluon has only two physical states. A po- 
larization vector on the other hand has four com- 
ponents but with two of the degrees of freedom 
implicitly removed by gauge invariance. It is this 
mismatch between the number of physical states 
and the number of components of the polariza- 
tion vectors that is at the origin of the apparent 
complexity. A much better representation of the 
amplitude may be found by making use of helicity 
or circular polarization. Using its most primitive 
implementation, one would replace each polariza- 
tion vector with an explicit circular polarization: 

£+ = (0,l,+z,0), = (0,l,-z,0). (9) 

However, this representation is not Lorentz co- 
variant. Much better Lorentz covariant represen- 
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Figure 10. The result of a brute force evaluation of the five-gluon tree amplitude on a computer in an 
unreadable font. The result shown here is expressed in terms of in terms of dot products of momentum 
and polarization vectors and is the coefficient of a single color factor. 



tations have been invented making use of 

spinor representations that are more convenient 
to use. This is sometimes referred to as 'Chinese 
Magic' in recognition of the especially beautiful 
formulation due to Xu, Zhang and Chang [3^ . 

In any case, if we make use of helicity, the com- 
plete physical content of Figure [l^ can instead be 
written much more neatly in terms of the color 
stripped helicity amplitudes: 



4+ 


5+) 


= 0, 


4+ 


5+) 






yS23S34S45S51 


4+ 


5+) 


_ ( 43 




\ S23S34S45S5I 



1/2 



(10) 



1/2 



where the = ipi+Pj)'^ are the kinematic invari- 
ants appearing in the ampUtude and the i and j 
are the labels distinguishing the different gluons. 
The + and — superscripts for each particle i la- 
bel the helicity (in a convention where all legs are 
outgoing). 

Given the remarkable cleanup of the ampli- 



tude, it is perhaps not surprising that helicity 
representations have allowed seemingly impossi- 
ble computations to proceed. Perhaps the most 
famous example of this is the n-point general- 
ization of eq. (10) for special helicity configura- 
tions, obtained by Parke and Taylor js^. They 
conjectured that in QCD the leading order color- 
stripped n-gluon amplitudes satisfy, 

A5(l±,2+,3+,4+,...,n+) = 0, 
A5(l-,2-,3+,4+,...,7i+) 



(11) 



*12 



1/2 



S23S34S45 • ■ • Snl 



, (12) 



for the special helicity configurations indicated by 
the plus and minus labels. Using standard Feyn- 
man rules one would need to evaluate an infinite 
number of diagrams to obtain this. However, us- 
ing helicity methods together with recursive tech- 
niques, a proof that the Parke- Taylor ampli- 
tude are correct was presented soon after the orig- 
inal conjecture. Helicity methods have also been 
used for constructing infinite sequences of one- 
loop amplitudes, again for cases of special helicity 
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configurations |^,^. As a more recent example, 
the first computation of a two-loop 2 — > 2 scatter- 
ing amplitude in QCD also relied on helicity 
methods. 



3.2. Applications of Supersymmetry 

Supersymmetry has a long history of use in 
constructing models of possible extensions to the 
Standard Model. Another application of super- 
symmetry is providing strong checks on non- 
trivial calculations in Standard Model physics. 
The high degree of symmetry renders the dy- 
namics of supersymmetric theories far simpler 
than those of non-supersymmetric ones used in 
describing the Standard Model. Calculations in 
supersymmetric theories can therefore serve as 
'toys' for devising methods for dealing with more 
complicated Standard Model calculations. More- 
over a QCD calculation is a close relative of a 
similar calculation in a supersymmetric version of 
QCD where the quarks are replaced by gluinos. 
This allows one to perform checks of QCD calcu- 
lations by performing minor alternations to the 
calculations so that they become supersymmet- 
ric. 

A particularly striking example of the use of su- 
persymmetry as a check is from the computation 
of the four-loop QCD /3-function performed by 
van Ritbergen, Vermaseren and Larin |^^. The 
(3 function is a fundamental quantity controlling 
the running of the coupling constant. This quan- 
tity was already known earlier in a supersymmet- 
ric version of QCD However, due to certain 
technical complications having to do with the dif- 
ferent ways that the theories were rendered finite, 
the two computations could not be compared di- 
rectly. Nevertheless a non-trivial cancellation of 
a particular color factor occurs when the QCD 
result is altered to be supersymmetric, in agree- 
ment with the expectation from the supersym- 
metric version. This suggests that the four-loop 
calculation of the QCD /3-function is indeed cor- 
rect. As a more recent example, supersymmetry 
identities have also been applied as checks on two- 
loop calculations of scattering amplitudes . 



3.3. Quantum Gravity Example of Growth 
of Difficulty 

It is not difficult to find problems in quantum 
field theory which cannot possibly be solved by 
brute force alone. In many cases systematic al- 
gorithms for evaluating Feynman diagrams are 
not known, but even where algorithms are known, 
problems can easily surpass the capacity of any 
computer. One example of a particularly inter- 
esting problem which possesses this property is 
the ultraviolet divergences of quantum gravity. 

The often-repeated statement that quantum 
mechanics and General Relativity are incompati- 
ble arises from the bad ultra-violet divergences of 
(super) gravity theories.^ There are strong argu- 
ments that this is in fact correct based on power- 
counting of the gravity Feynman diagrams (see, 
e.g., ref. @). 

For the case of pure Einstein gravity, these 
arguments were confirmed by an explicit two- 
loop computation performed by Sagnotti and Go- 
roff |47| . However, no proof exists for supergrav- 
ity theories. The first potential divergence al- 
lowed by supergravity occurs at three quantum 
loops. One of the Feynman diagrams which would 
need to be evaluated to obtain the divergence is 
shown in Figure ^ A naive estimate of the total 
number of terms in this diagram, based on the 
structure of the propagators and vertices yields 
approximately 10^^ terms. If one could evalu- 
ate terms at the rate of a billion per second, it 
would take approximately 20,000 years to extract 
the ultra-violet divergence of just this single dia- 
gram. It is no surprise that this computation has 
not been attempted. 

There are, however, methods which can bring 
this problem into the realm of possible calcu- 
lations. A trick is to map the gravity calcula- 
tion in products of gauge theory calculations us- 
ing certain string theory relations between am- 
plitudes |4^,^. The net effect of this is that an 
impossible problem can be mapped into an ex- 
tremely difhcult problem, which also has not been 
solved, as yet. This technique has, however, al- 

^ There actually is no incompatibility in perturbation the- 
ory, up to an energy scale of 10^^ GeV, if General Rela- 
tivity is viewed as an effective field theory jZs) instead of 
as a fundamental theory. 



8 



Figure 11. A sample diagram whose divergence 
part would need to be evaluated in order to deter- 
mine the ultra-violet divergence of a supergravity 
theory. The lines represent graviton propagators 
and the vertices three-graviton interactions. 

ready been used to show that at least for the case 
of maximally supersymmetric gravity the onset of 
divergences is delayed until at least five quantum 
loops 

4. STATUS OF LOOP CALCULATIONS 

Before surveying the main advance since the 
last ICHEP conference, it is useful to survey the 
status of quantum loop calculations. Here we do 
not discuss tree-level calculations which have also 
seen considerable progress over the years. 

4.1. Status of one-loop calculations 

In 1948 Schwinger dealt with one-loop three- 
point calculations such as that of the anoma- 
lous magnetic moment of leptons described in 
Section ||. It did not take very long be- 
fore Karplus and Neuman calculated light-by- 
light scattering in QED in their seminal 1951 
paper [^. In 1979 Passarino and Veltman pre- 
sented the first of many systematic algorithms for 
dealing with one-loop calculations with up to four 
external particles, leading to an entire subfield de- 
voted to such calculations. Due to the complexity 
of non-abelian gauge theories, however, it was not 
until 1986 that the first purely QCD calculation 
involving four external partons was carried out in 
the work of Ellis and Sexton . 

The first one-loop five-particle scattering am- 
plitude was then calculated in 1993 by Lance 
Dixon, David Kosower and myself ||5^ for the 
case of five-gluon scattering in QCD. This was 
followed by calculations of the other five-point 
QCD subprocesses Q, with the associated phys- 



ical predictions of three-jet events at hadron col- 
liders appearing somewhat later (s^,^ . A num- 
ber of other five-point calculations have also been 
completed. One example of a state-of-the-art five- 
point calculation was presented in a parallel ses- 
sion by Doreen Wackeroth , who described the 
calculation of pp ttH at next-to-leading order 
in QCD ||5^. This process is a useful mode for 
discovering the Higgs boson as well as measure- 
ment of its properties. Other examples are NLO 
calculations for e+e^ 4 jets [ ^9|j60| , |6H , Higgs 
-I- 2 jets 16^, and vector boson -f 2 jet produc- 
tion [ p9|j63[ |, which is also important as a back- 
ground to the Tevatron Higgs search, if the jets 
are tagged as coming from b quarks. 

Beyond five-external particles, the only calcu- 
lations have been in special cases. By making 
use of advanced methods, for special helicity con- 
figurations of the particles, infinite sequences of 
one-loop amplitudes with an arbitrary number 
of external particles but special helicity configu- 
rations have been obtained in a variety of the- 
ories |39| , ^ . For the special case of maximal 
supersymmetry, six-gluon scattering amplitudes 
have been obtained for all helicities |^^. There 
has also been a recent calculation of a six-point 
amplitude in the Yukawa model |Q , as well as re- 
cent papers describing properties of six-point in- 
tegrals Q . These examples suggest that that the 
technical know-how for computing general six- 
point amplitudes is available, though it may be 
a rather formidable task to carry it through. An 
efficient computer program for dealing with up to 
three jets at hadron colliders now exists [Q, sug- 
gesting that it would be possible add one more 
jet, once the relevant scattering amplitudes are 
calculated. This would then give a much bet- 
ter theoretical handle on multi-jet production at 
hadron colliders. 

4.2. Status of Higher Loop Computations 

Over the years, an intensive effort has gone 
into calculating higher loop Feynman diagrams. 
A few samples of some impressive multi-loop cal- 
culations are: 

• The anomalous magnetic moment of lep- 
tons, already described in Section ||. 
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hadrons) /cr(e+e 



which has been calculated through 
0{a^g) This quantity is one of the cor- 
nerstones demonstrating that QCD is the 
correct theory of strong interactions. 

• The four-loop QCD /3 function computed by 
van Ritbergen, Vermaseren and Larin p2[ . 
In the course of this computation approxi- 
mately 50,000 Feynman diagrams were eval- 
uated, of which a single one is displayed in 
Figure p. 




Figure 12. One of the 50,000 Feynman diagrams 
appearing in the computation of the four-loop 
QCD 13 function. 



These types of precision calculations have been 
crucial for comparing the Standard Model to ex- 
periment |67j. However, except for very special 
cases, all two and higher-loop Feynman diagram 
computations prior to a few years ago contained 
either one or no kinematic variables. 

5. ADVANCES OF PAST TWO YEARS 

The key advance in perturbative computations 
over the last two years is our ability to perform 
two-loop computations with more than a single 
kinematic variable. Previously, the only such cal- 
culations were for special cases of maximal super- 
symmetry [Q|9|. The pace of progress, due to a 
large extend to the influx of energetic young peo- 
ple into the field, is such that even though two 
years ago it was not known how to do such calcu- 
lations for more general processes, now they are 
commonplace. A sample diagram illustrating the 
types of processes that can now be computed is 
depicted in Figure ^ 

Every step in the construction of a physical 
cross-section involving a two-loop amplitude has 



o 
o 
o 



AAAAA/ 



Figure 13. An sample process that has been eval- 
uated. The curly lines represent gluons, the wavy 
ones photons and the straight ones quarks. 



serious technical challenges. Although some of 
these difficulties have not been completely re- 
solved, all aspects of the problem have seen re- 
markable progress over the last few years. 

5.1. Loop Integrals 

A crucial ingredient in this progress has been 
the breakthrough in obtaining the integrals 

p massless 2^2 
2| . For example, 



needed for computing two-loo 
scattering amplitudes 1 69 , 7| 



consider the scalar double box integral. 



1 



1 

{p- ki - {q - ki)'^ {q ~ ks - k^Y 



(13) 



which may be represented by the left diagram in 
Figure 14. This scalar integral was a non-trivial 
challenge, until it was evaluated by Smirnov |69| 
in terms of standard polylogarithms. Shortly 
thereafter, the non-planar scalar integral appear- 
ing on the right-hand-side of Figure |lj as well 
as other relevant integrals were evaluated ||7^,0 . 
Moreover, powerful new methods for reducing 
general tensor integrals into a basis of known inte- 
grals were developed [Q, making use of integra- 
tion by parts j73| and Lorentz invariance |74|. A 
number of important further improvements have 
also been made allowing for more systematic eval- 
uations |7^. The final expressions obtained with 
these methods are all in terms of standard func- 
tion such as logarithms, polylogarithms as well 
as generalizations of these functions that are also 
amenable to efficient numerical evaluation. 
5.2. Amplitude Calculations 

The first Standard Model calculation of a 2 ^ 2 
scattering amplitude was that of gluon scattering 
in QCD for a special hehcity configuration pl| . 
By making use of the loop integration break- 



10 




Figure 14. Planar and non-planar double box di- 
agrams. 



through, since then many new two-loop ampli- 
tudes appeared: 

• Bhabha scattering in QED, e+e^ e+e^, 
in the ultra- relativistic limit ||lj]. Bhabha 
scattering is useful for monitoring luminos- 
ity at an e'^e~ collider. 

• All two-loop 2 — > 2 QCD processes ||l^, 
[Tl| . These are essential ingredients for con- 
structing NNLO jet programs for hadron 
colliders. 

• Light-by-light scattering, 77 77, in the 
ultra-relativistic limit [ j76[ , which is of inter- 
est for future photon-photon colliders. 

• Gluon fusion into a photon pair, gg 

77 Q . As described in Section ^, this sub- 
process is important as a background to 
Higgs production when a light Higgs decays 
into a pair of photons. 

• Two loop QED and QCD corrections to 
massless fermion boson scattering f77| : 
qq — > 77, qq g7, and e+e" 77. 

• e+e^ — > 3 partons |l^. This is a key 
process for a future high energy e^e~ col- 
lider and it will allow precision measure- 
ments of strong coupling processes at the 
1% level |7|. 

• Deeply inelastic scattering 2jet production 
and pp^W,Z + 1 jet |l|,|9[. 

As more integrals including ones with mas- 
sive legs are worked out |80[| , other amplitudes in 
QED, electroweak theory, heavy quark physics, 
etc. will surely also be evaluated. Such calcula- 
tions are just the initial steps since they are for 



matrix elements and not for physical predictions. 
As yet only the 77 — > 77 and gg — > 77 ampli- 
tudes have been implemented in phenomenologi- 
cal studies [f76|,p|. 

5.3. Infrared Divergences 

The source of much grief in perturbative cal- 
culations are infrared divergences. At the end, 
for any appropriately defined physical quantity 
these divergences must cancel |8^. However, at 
intermediate steps of calculations, before the var- 
ious contributions are combined to form a phys- 
ically meaningful quantity, severe divergences 
arise from soft (i.e. pi 0) or coUinear particles. 
At a given order or perturbation theory, a phys- 
ical quantity typically involves both virtual and 
real emission contributions that are separately in- 
frared divergent. 




fa) (h) (c) 



Figure 15. The various contributions at NNLO 
to the squared matrix elements: (a) virtual, (b) 
single real emission and (c) double real emission. 



At two-loops it becomes much more difficult 
to deal with these divergences due to their intri- 
cate structure. An important result is Catani's 
Magic Formula' which specifies the precise 
form of the two-loop divergences. The 'magic' 
is that it was obtained before any explicit calcu- 
lations were possible. Feynman supposedly once 
remarked that before one calculates one should 
know the result. Needless to say, this ideal is dif- 
ficult to achieve, but Catani's Magic Formula is 
a partial implementation in the sense that the in- 
frared divergences that must emerge are known 
prior to beginning a calculation. Catani's Magic 
Formula then provides important guidance in any 
explicit two loop calculation. It also provides a 
good way for organizing amplitudes into diver- 
gent parts which ultimately drop out of physical 
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quantities and the finite contributions which de- 
termine the physics. 

In order for the advances in computing scat- 
tering amphtudes to be used in comparisons of 
theory to experiments, the various contributions 
need to be combined into a numerical program for 
producing the plots comparing theory to experi- 
ment. In general, this is non-trivial because of the 
infrared divergences plaguing intermediate steps 
of the calculations. At NLO general methods for 
dealing with the infrared singularities exist Q, 
but at NNLO as yet there is no complete solution. 
There is, however, a partial understanding. In 
particular, the universal behavior of amplitudes 
entering into NNLO calculations as the momenta 
of partons become either soft or coUinear has been 
worked out js^,^. Moreover, in certain special 
cases, such as Drell-Yan p6[, inclusive Higgs pro- 
duction at hadron colliders |7|,^ and e^e~ 7-I-I 
jet this has been worked out. It is now clear 
that more general methods for dealing with the 
infrared divergences at NNLO will be constructed 
in the near future. 



5.5. Problems with multiple mass scales. 

For problems with multiple mass scales, in re- 
cent years there has also been quite a bit of 
progress. A powerful general strategy is known 
as 'the strategy of regions' |^,p6 33 . The basic 



trick makes use of dimensional regularization al- 
lowing one to series expand expressions in various 
kinematic limits or regions without the need for 
keeping track of the boundaries between regions. 
The key progress is that now there is a universal 
method for dealing with this, instead of case-by- 
case analysis depending on boundaries between 
the regions. This has been applied to a large 
number of problems | |9^ , ranging from atomic 
physics, heavy quark production and decay, large 
electroweak logarithms and very recently lattice 
gauge theory. As one example from this con- 
ference, the electroweak logarithms described by 
Kiihn were obtained using the strategy of re- 
gions. 



5.4. Parton Distribution Functions 

In order to have true NNLO predictions for 
observables at hadron colliders a crucial ingredi- 
ent is NNLO parton distribution fimctions. At 
NLO the evolution kernels needed for the DGLAP 
equations Q were calculated long ago Q . The 
calculation of the next order has been a diffi- 
cult challenge since then, but recent years have 
seen some impressive progress. The problem is 
partially solved pO| , and very good approximate 
solutions based on computing Mellin moments 
have also been obtained These are starting 
to be implemented in global fits by the MRST 
group There is also some related work on 

NNLO quark and gluon distributions inside pho- 
tons instead of protons. We can anticipate 
that it will not be long before all standard par- 
ton distribution functions are fully implemented 
at NNLO. (A description of other important par- 
ton distribution function issues, such as uncer- 
tainties, may be found in the talks of Giele and 
Stump [|4|.) 



5.6. Numerical Methods 

Another area where there has been recent 
developments is in purely numerical techniques 
for evaluating Feynman diagrams. An impor- 
tant recent advance in numerical methods was 



made by Ghinculov and Yao |99 , who produced 
a method for dealing with two-loop problems 
where the infrared divergences are mild, as typi- 
cal in problems in electroweak theory. An exten- 
sion of this was described in the parallel session 
talk of Passarino |100[ , targeting state-of-the-art 
electroweak calculations. In QCD where there 
are severe infrared divergence, anothe r nu meri- 
cal method, by Binoth and Heinrich |101|, has 
been applied as an important check on state-of- 
the-art analytical two-loop calculations. At this 
conference Dave Soper described a purely numer- 
ical method for massless QCD amplitudes [102|. 
So far, this latter method has been used to re- 
produce the seminal result for e^e~ — > 3 jets at 
NLO obtained by Ellis, Ross and Terrano |103| in 
1981, but it is not yet clear if this method will be 
applied to state-of-the-art calculations. 
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6. SAMPLE APPLICATIONS OF THE 
BREAKTHROUGH 

In the month prior to the conference, the first 
two applications of the new advances in two-loop 
calculation to specific problems of interest for 
collider physics appeared |^,^. These first ex- 
amples focus on Higgs physics in part because 
of its importance for the collider physics pro- 
gram, but also because the infrared singularities 
for these cases are relatively straightforward to 
handle. Once the general problem of infrared 
divergent phase space integrals is solved, many 
more applications will appear. 

The first of these applications provides an im- 
proved understanding of the QCD background to 
Higgs production and decay for the case of a light 
Higgs. The second application, which appeared 
just two weeks prior to the conference, provides 
the exact result for inclusive Higgs production ||^ 
at NNLO. This result is in perfect agreement with 
an earlier series expansion due to Harlander and 
Kilgore [Q and in very good agreement with an 
even earlier approximate solution [104 105||. 



6.1. Inclusive Higgs Production 

In his talk, Frixione ]l3| described the impor- 
tant quest to obtain reliable predictions for Higgs 
production at the Tevatron and LHC. Previously, 
inclusive Higgs production had been calculated 
through NLO |106], leaving a substantial theoret- 
ical uncertainty in the result. The new NNLO cal- 
culations significantly reduce the theoretical un- 
certainties. 

Because of the large Yukawa coupling of the 
Higgs to top quarks and the fact that the LHC 
acts as a glue factory the dominant Higgs pro- 
duction mechanism LHC is though gluon fusion 
via a top loop. For a light mass Higgs boson, 
preferred by precision electroweak measurements, 
the top loo p ca n me replaced by a simpler effec- 
tive vertex [107|, as depicted in Figure |l^. In this 
case, the effective vertex is an excellent approxi- 
mation Ip^ . 

Using the effective Higgs vertex, the two-loop 
virtual contributions of the type in Figure p^(a) 
were worked out recently by Harlander |109|. The 
especially difficult part of the calculation is deal- 




Figure 16. For a light Higgs the top quark loop 
can be replaced by an effective vertex. The curly 
lines represent gluons, the solid ones the top 
quark, and the dashed ones the Higgs boson. 



ing with the IR singular integration over the dou- 
ble real emission phase space encountered when 
evaluating contributions arising from diagrams of 
the type shown in Figure |l^(c) . This was solved 
recently by Harlander and Kilgore who con- 
structed a series expansion of the NNLO inclusive 
cross section. 

The next advance was by Anastasiou and Mel- 
nikov who applied the new loop integration tech- 
nology outlined in Section 5.1 to obtain the exact 
phase space integration over the various contribu- 
tions including the double real emission ones 
The essential trick for doing this is based on 
unitarity, as encoded in the Cutkosky cutting 
rules I lie]. The relationship between phase space 
and Feynman loop integrals is a rather useful 
trick, having, for example, also been used to 
obtain a number of two-loop scattering ampli- 
tudes 0||jl^. With this trick it is possible to re- 
place the phase space integration, including those 
for double real emission, with standard Feynman 
loop integrals. Once the problem is expressed in 
terms of standard Feynman integrals the new in- 
tegration methods can be applied directly. So 
far this trick has been applied to totally inclu- 
sive cross-sections, but it is very likely that it can 
be extended to at least certain differential cross- 
sections 
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The basic result of both these calculations is 
summarized in Figure |l^, obtained from ref. 
The perturbative expansion clearly stabilizes at 
NNLO and is reliable, given the modest shift be- 
tween the NLO and NNLO results. Further re- 
finements to the NNLO results should also be 



forthcoming [105|. 
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Figure 17. Sample diagrams representing (a) virtual, (b) single emission and (c) double emission contri- 
butions. The curly lines represent gluons, the dashed ones Higgs bosons, and the solid circle the effective 
vertex coupling gluons to the Higgs. 



,a(pp^H+X)[pb] 



Vs"= 14 TeV 
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Figure 18. The LO, NLO and NNLO predictions 
for inclusive Higgs productions at the LHC as a 
functions of the Higgs mass (with factorization 
and renormalization scales set equation to the 
Higgs mass). 



6.2. The Di-Photon Background to Higgs 
Decay 

The second application of the breakthrough 
that we describe is for the background to ob- 
serving a light {Mh < 140 GeV) Higgs boson 
at the LHC. For a light Higgs, a key search mode 
at the LHC involves Higgs production via gluon 
fusion followed by the rare decay into a pair of 
photons. (For a recent survey of ways to detect 
a light Higgs boson at the LHC see, for exam- 
ple, ref. |ll2[| .) Although the branching ratio 
is tiny, this mode is relatively clean due to the 
excellent mass resolution of the LHC detectors. 
This allows the background to be measured ex- 
perimentally and subtracted from a putative sig- 



nal (TT3|jTT|,|TT§. Nevertheless, it is stiU impor- 
tant to have robust theoretical predictions in or- 
der to systematically study the dependence of the 
signal relative to the background so as to optimize 
Higgs search strategies. Given that for the case 
of a light Higgs boson it will take about two years 
of running at the LHC before the Higgs signal is 
pulled out of the background, there is good moti- 
vation for wanting to optimize search strategies. 

The background is composed of two pieces. 
The 'reducible' background is where photons are 
faked by jets, or more generally by hadrons, es- 
pecially 7r"s. This background can be suppressed 
efficiently by photon isolation cuts, where events 
are rejected based on the amount of hadronic en- 
ergy near the photons. The second source of back- 
ground is from the underlying QCD process where 
quarks emit photons either directly or through 
fragmentation. The process pp JjX proceeds 
at lowest order via the quark annihilation subpro- 
cess qq —^ 77, which is independent of the strong 
coupling a^. One of the Feynman diagrams de- 
scribing this process is shown in Figure |l^ The 
next-to-leading-order (NLO) corrections to this 
subprocess have been incorporated into a number 
of Monte Carlo pro gram s |116|| , the most up-to- 
date being DIPHOX JUTt . 



,1 



Figure 19. A leading order diagram contributing 
to di-photon production. 
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Although the gluon fusion contribution de- 
picted in Figure ^ is formally of higher or- 
der in the QCD coupling, it is enormously en- 
hanced by the fact that the distribution in the 
proton becomes very large at small x. This 
makes formally higher order corrections involv- 
ing gluon initial states very significant for the 
production of low-mass systems (< 200 GeV) at 
the LHC. The net result is that the gluon fu- 
sion contribution to pp — > ^^X is comparable 
to the leading-order quark annihilation contribu- 



tion [116117,118,1191 



5'^ 



9 
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Figure 20. A leading order diagram contributing 
to the gluon fusion into two photons. 



loop contributions of the type in Figure |2l|(a) 
must be combined with the real emission contri- 
butions in Figure |l|(b) in such a way that all 
infrared singularities cancel. In this case the in- 
frared singularities are all in a form which can be 
handled using standard NLO methods . Then 
these new contributions must be combined with 



the previously obtained ones |117|. 

A sample result taken from ref. 1^ is shown in 
Figure |2^. The conclusion from this plot is that 
in the region of interest the NLO corrections to 
the gg — > 77 subprocess have a relatively mod- 
est effect on the total irreducible di-photon back- 
ground to the Higgs search. Indeed the K fac- 
tor (ratio of NLO to LO contributions) for this 
subprocess is significantly smaller than previous 
estimates somewhat enhancing the statistical sig- 
nificance of the Higgs signal compared to previ- 
ous estimates. Moreover, the relatively small shift 
from LO to NLO suggests that the subprocess is 
under adequate theoretical control |9| . 



To reduce the uncertainty on the total 77 pro- 
duction rate, one therefore needs to calculate the 
gg 77 subprocess at its next-to-leading-order, 
even though it is formally N'^LO as far as the 
whole process pp — > 77X is concerned. This in- 
volves the diagram of the type shown in Figure |2l|. 
The two-loop virtual contributions (a) were com- 
puted recently in ref. while the one- loop real 
emission contri buti on (b) is obtained from a per- 
mutation sum [120 1 over contributions to the one- 
loop five-gluon amplitude |53] . 
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Figure 21. Sample NLO diagrams contributing 
to gluon fusion into two photons: (a) virtual and 
(b) real emission contributions. 



a 

•a 



b 
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pp ~* 77X Cross Section at LHC 
I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' 
NLO 




Figure 22. The pp — > 77 AT cross-section as a func- 
tion of invariant mass of the di-photon pair at the 
LHC. The lower curve is obtained using DIPHOX. 
The upper curve is the result of adding in the 
extra NLO contributions to gluon fusion. 



To obtain the physical cross-section the two- 



A crucial question is whether we can improve 
the situation by finding appropriate search strate- 
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gies which enhance the signal over the back- 
ground. Some initial studies of this question at 
the parton level were performed in ref. but 
more realistic studies are required, including de- 
tector effects illSll and the reducible 7r° back- 



ground contributions |114,115 121] 



7. SUMMARY AND OUTLOOK 

A very recent example illustrating the impor- 
tance of precision perturbative calculations is the 
impressive measurement of the anomalous mag- 
netic moment of the muon by the Brookhaven 
g~2 collaboration. The first part of this talk out- 
lined the intense theoretical effort that has gone 
into obtaining the Standard Model prediction at 
a precision matching the one of the experiment. 

This talk also highlighted the recent rapid pace 
of progress in our ability to compute quantities of 
interest in perturbative quantum field theory. In 
particular, two-loop calculations involving up to 
four external particles can now be performed in 
rather general cases. This breakthrough is the 
culmination of years of effort. Parallel advances 
in calculations of the DGLAP evolution of parton 
distribution functions have also been made. 

So far the breakthrough had been applied 
mainly to QCD and to massless QED, but as more 
types of loop integrals with masses are calculated 
one can expect applications also to electroweak 
theory, heavy quark physics, QED with massive 
particles, and so forth. Two very recent phe- 
nomenological applications of the advance are the 
calculation of the exact next-to-next-to-leading 
order inclusive Higgs production cross-section at 
hadron colliders and the calculation of the QCD 
background to Higgs production and decay into 
a photon pair at the LHC, for the case of a 
light Higgs. Once general algorithms are set up 
for dealing with infrared divergent phase space 
at next-to-next-leading order, many more phe- 
nomenological studies will appear. 

Given the influx into the field of the many 
talented energetic young people who have con- 
tributed greatly to the advances described here, 
we can be optimistic that the rapid pace of 
progress will continue for the foreseeable future. 
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